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Abstract 

Traveling wave solutions of degenerate coupled UKdV equations are studied. Due 
to symmetry reduction these equations reduce to one ODE, (/')^ = Pn{f) where Pn{f) 
is a polynomial function of / of degree n = £ + 2, where £ > 3 in this work. Here I is the 
number of coupled fields. There is no known method to solve such ordinary differential 
equations when i > 3. For this purpose, we introduce two different type of methods 
to solve the reduced equation and apply these methods to degenerate three-coupled 
KdV equation. One of the methods uses the Chebyshev’s Theorem. In this case we find 
several solutions some of which may correspond to solitary waves. The second method 
is a kind of factorizing the polynomial Pn{f) as a product of lower degree polynomials. 
Each part of this product is assumed to satisfy different ODEs. 

Keywords: Traveling wave solntion. Degenerate KdV system, Chebyshev’s theo¬ 
rem, Alternative method. 


1 Introduction 

The system of degenerate coupled multi-field KdV equations is given as PP-i] 

3 2 

Uf ^UUx T Qx 

2 2 1 ^ 2,3 

Qt = q Ux+ ^UQx + Qx 

f_-\ tf_i 1 /j_-| 

Qt = q Ux + -UQx + Vx 
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where = u and = v. In n-UI. it was shown that this system is also a degenerate KdV 
system of rank one. In a previous work [S] we focused on the equation (1.1) for i = 2. We 
reduced this equation into an ODE (/')^ = P^i^f) where P^if) is a polynomial function of 
degree four. We analyzed all possible cases about the zeros of P^if). Due to this analysis 
we determined the cases when the solution is periodic or solitary. When the polynomial has 
one double /2 and two simple zeros /i, /s with /i < /2 < /s, or one triple and one simple 
zeros, the solution is solitary. Other cases give periodic or non-real solutions. By using the 
Jacobi elliptic functions [9], we obtained periodic solutions and all solitary wave solutions 
which rapidly decay to some constants, explicitly. We have also shown that there are no real 
asymptotically vanishing traveling wave solutions for i = 2. Indeed we have the following 
theorem for the degenerate coupled £-KdV equation, i > 2. The degenerate coupled £-KdV 
equation can be reduced to the equation 


{fr = ( 1 . 2 ) 

by taking i functions as u{x,t) = q^{x,t) = /(O, = /2(0, • • •, v{x,t) = q^ = /^(O, 

where ^ = x — ct in (1.1). Here P£+ 2 (/) is a polynomial of / of degree i + 2. If we apply the 
asymptotically vanishing boundary conditions to (1.2), we have 


(fr = Bf‘(f + 2 c)', 


(1.3) 


Theorem 1.1 When i = odd, we have P > 0, and then the degenerate coupled i-KdV 
equation has real traveling wave solution with asymptotically vanishing boundary conditions, 
but when i = even, the constant B < 0. Hence the equation does not have a real traveling 
wave solution with asymptotically vanishing boundary conditions. 


Proof. Consider the degenerate coupled PKdV equation (1.1). Let u{x,t) = q^{x,t) = f{^), 
g^(x, t) = /2(0) • • •) t) = = fiiO where ^ = x — ct. By using the first £ — 1 equations 

we obtain all functions i = 2, 3,..., £ as a polynomial of /(O- We get 


m) = Q2(f) = -cf-^f + d, = -C,2f + A,(f) 

— Q:i{f) — B (c^ — di)f + d2 

= + ^2{f) 

fiiO = Qiif) = ~ 2 ^-^^ ~ ^ ~ 

= + As{f) 


feiO = Qe{f) = {-lY^Wf + Ae.^if), 
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where aj > 0 are constants, Ai{f) and Qj{f) are polynomials of / of degree i and j, 
i = 1, 2, 1, j = 2,3,..., i, respectively. 

Now use the £th equation. We obtain 





l,8Ai-Yf) 8A,.,(f) 


dx 


dt 


(1.4) 


Integrating above equation once we get 


1 

4 


/" 




By using /' as an integrating factor, we integrate once more. Finally, we obtain 


U'f = (- 1 ) 


£-1 




(£+!)(£+ 2 ) 


+ ^) + R^+iU) = + Re+i{f), (1.5) 


where Ri{f) is a polynomial of / of degree i, i = i,£ + 1. When i is odd, the coefficient 
of that is B, is positive, and when i is even, B is negative. Applying asymptotically 
vanishing boundary conditions, we get (/O^ = Bp{f + 2cY, where sign(i?) = (—1)^“^. Hence 
for £ = odd, the degenerate coupled £-KdV equation has real traveling wave solution with 
asymptotically vanishing boundary conditions, but when £ = even, it does not. □ 


In this work, we study the equation (1.1) for £ 


3 which is 


3 

Uf 2^^^: ffi Vx 

1 

Vt = VUx + -UVx + OJx 

1 1 

OJl ^^xxx ffi ^Ux ~U^xi 


( 1 . 6 ) 

(1.7) 

( 1 . 8 ) 


in detail. It is clear from Theorem ] 1.1 that unlike the case £ = 2, 
solution with asymptotically vanishing boundary conditions in £ 


we have real traveling wave 
= 3 case. 


Here the system (1.6)-(1.8) reduces to a polynomial of degree hve. 


{f'Y ~ ^ + (6c^ — 2di)p + 4(c^ — cdi + d2)P + Sd^f + 8d4 — PYf)i (1-9) 


where fp) = u{x,t), and c,di,d 2 ,d^,d 4 are constants. When the degree of the polynomial 
in the reduced equation is equal to hve or greater it is almost impossible to solve them. As 
far as we know there is no known method to solve these equations. We shall introduce two 
methods to solve such equations. The hrst one is based on the Chebyshev’s Theorem [10] 
which is used recently to solve the Einstein held equations for a cosmological model HH, 
[T2|. By using the Chebyshev’s theorem we give several solutions of the reduced equations for 
£ = 3 and also for arbitrary £. The second method is based on the factorizing the polynomial 
F£+ 2 (/) as product of lower degree polynomials. In this way we make use of the reduced 
equations of lower degrees. We have given all possible such solutions for £ = 3. 
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The layout of our paper is as follows: In Sec. II, we study the behavior of the solutions in 
the neighborhood of the zeros of P 5 {f) and discuss the cases giving solitary wave solutions. 
In Sec. Ill, we hnd exact solutions of the system (1.6)-(1.8) by a new method proposed for 
any £> 3 which uses the Chebyshev’s Theorem and analyze the cases in which we may 
have solitary wave and kink-type solutions. In Sec. IV, we present an alternative method. 
Particularly, we hnd the solutions of the system (1.6)-(1.8). Here we obtain many solutions 
for £ = 3 including solitary wave, kink-type, periodic, and unbounded solutions. We present 
some of them in the text but the rest of the solutions are given in the Appendices A and B. 


2 General Waves of Permanent Form for (i = 3 ) 


2.1 Zeros of P 5 {f) and Types of Solutions 


Here we will analyze the zeros of P^{f) in (1.9) 


(i) If /i = /(^i) is a simple zero of P^^f) we have P^{fi) = 0. Taylor expansion of P^^f) 
about /i gives 


iff = P^{h) + P',{h){f - h) + 0{{f - hf) 
= P',{h)U-h) + o{{f-hf). 


From here we get = 0 and f\ii) = P^{fi)/2. Hence we can write the function /(^) as 

/K) = /Ki) + K-?.)/'(?.)+ ^K-?.)V"Ki) + o((f-6)t 

= h + + ( 2 . 1 ) 

Thus, in the neighborhood of .^ = .^i, the function f{^) has local minimum or maximum as 
P^{fi) is positive or negative respectively since = P'^{f\)/2. 


(ii) If fi = /(^i) is a double zero of P^^f) we have P 5 {fi) = P^ifi) = 0. Taylor expansion 
of P^if) about fi gives 


iff = n(/i) + i"5(/i)(/-/i) + o(/-/irWi) + o((/-/in 




( 2 . 2 ) 


To have real solution /, we should have P^{fi) > 0. From the equality (2.2) we get 
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which gives 


(2.3) 

where a is a constant. Hence / —)■ /i as =Fcx). The solntion / can have only one peak 
and the wave extends from — cx) to cx). 


(iii) If /i = /(^i) is a triple zero of P^if) we have P 5 (/i) = P^ifi) = = 0- Taylor 

expansion of P 5 {f) abont fi gives 

(ff = Wi) + f?(/i)(/-/i) + i(/-/i)"P5(/i) + ^(/-/i)’ + 0((/-/i)-‘) 

= ^(/-/itWO + O))/-/!)-). (2,4) 

This is valid only if both signs of (/ — /i)^ and P^"{fi) are same. Hence, to obtain real 
solntion / we have the following two possibilities: 

1 ) (/-A)>0andPr(/i)>0, 

2 ) (/-A)<0andPr(/i)<0. 


If (/ — /i) > 0 and P^'ifi) > 0 then we have 


which gives 


f ^ ±—(f - 


/(?) ~ /i + 


± ^ vPU)? + «i) 

where ai is a constant. Thus / —)■ /i as ±cx. 

Let (/ — /i) < 0 and P^”{fi) < 0. Then 


2 ’ 


which yields 




/(O ~ /i - 


( + ^ 2 ) 

where 0:2 is a constant. Thus / —)■ /i as ±cx. 


2 ’ 


(2.5) 


( 2 . 6 ) 


(iv) If /i = /(&) is a quadruple zero of Ps)/), then we have Ps{fi) = P^(fi) = P('(/i) = 
P("(/i) — 0, In this ease Taylor expansion of Pe(f ) about /, gives 

(ff = ^(/ - /O-Pf)/.) + 0((/ - /,)=). (2.7) 
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This is valid only if > 0. Then we have 




which gives 


/(O ~ /i 


±^Vn^"(/i)^ + 7i 


iVeV • 

where 71 is a constant. Thus / —)■ /i as ±cxd. 

(v) If /i = /(^i) is a zero of multiplicity 5 of P^^f), then we have P^{f) 
is valid only if / — /i > 0. So we obtain the solution / as 

.4x1/3 1 


/4\ - 

/ = /l + (g) 


(±i + m,) 

where mi is a constant. Hence / —>■ /i as ± 00 . 


\ 2/3’ 


( 2 . 8 ) 

if-fif/2. This 
(2.9) 


2.2 All Possible Cases Giving Solitary Wave Solutions 

We analyze all possible cases about the zeros of P 5 {f) that may give solitary wave solutions. 
Here in each cases we will present the sketches of the graphs of P^if)- Real solutions of 
(y ')2 = P^^f) > 0 occur in the shaded regions. 

(1) One double and three simple zeros. 



Figure 1: Graphs of P^if) having one double and three simple zeros 

In Figure l.(b), /i, /a, and are simple zeros and /2 is a double zero. The real solution 
occurs when / stays between fi and /2 or /2 and /a. At /i, P^ifi) = > 0 hence graph 

of the function / is concave up at ^ 1 . At double zero / 2 , / —)■ /2 as ^ ± 00 . Hence we have 

a solitary wave solution with amplitude /i — /2 < 0. Similarly at /a, P^ifs) = < 0, 

hence graph of the function / is concave down at .^a- Therefore, we also have a solitary wave 
solution with amplitude 73-/2 >0. 

Now consider the graph (d) in Figure 1. For /a we have P^ifs) = > 0 thus graph of 

the function is concave up at .^ 3 . At double zero / 4 , / —)■ /4 as ±cxd. Hence we have a 
solitary wave solution with amplitude /a — /4 < 0. In other cases we have periodic solutions. 
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(2) Two double and one simple zeros. 




Figure 2: Graphs of P 5 {f) having two double and one simple zeros 

In Figure 2.(f), fi and /s are double zeros and /2 is a simple zero. The real solution occurs 
when / stays between /2 and /s. For /2 we have P^{f 2 ) = > 0 thus graph of the 

function is concave up at ^2- At double zero /s, / —>■ /s as ^ ±oo. Hence we have a 
solitary wave solution with amplitude /2 — /s < 0. 

Now consider the graph (g) in Figure 2. Here /2 and /s are double zeros and /i is a simple 
zero. The real solution occurs when / stays between /i and /2 or /2 and /s. For fi we have 
> 0 thus graph of the function is concave up at .^i. At double zero /2, / —)■ /2 
as ±oo. Hence we have a solitary wave solution with amplitude /i — /2 < 0. The other 
cases give kink, anti-kink type or unbounded solutions. 

(3) One triple and two simple zeros. 



Figure 3: Graphs of Pb{f) having one triple and two simple zeros 

Consider the graph (h) in Figure 3. Here /2 and /s are simple zeros and /i is a triple zero. 
The real solution occurs when / stays between /i and /2. For /2 we have P^{f 2 ) = < 0 

so graph of the function / is concave down at ^2- At triple zero /i, / —)■ /i as ^ ±oo. 

Hence we may have a solitary wave solution with amplitude /2 — /i >0. 

In the graph in Figure 3.(i), /i and /s are simple zeros and /2 is a triple zero. The real 

solution occurs when / stays between fi and /2. For fi we have P^{fi) = > 0 so 

graph of the function / is concave up at .^i. At triple zero /2, / —t /2 as ^ ±cxd. Hence we 
may have a solitary wave solution with amplitude /i — /2 < 0. The other case gives periodic 
solution. 
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(4) One quadruple and one simple zeros. 






Figure 4: Graphs of having one quadruple and one simple zeros 

In the graph Figure 4.(1), fi is a simple zero and /2 is a quadruple zero. The real solution 
occurs between /i and f^- At /i we have P^{fi) = > 0, so graph of the function / is 

concave up at ^i. At quadruple zero /2, / —)■ /2 as ^ ±cx). Hence we may have a solitary 

wave solution with amplitude /i — /2 < 0. The other case gives unbounded solution. 

(5) One double and one simple zeros. 

./l ,/2 

(n) 



Figure 5: Graphs of P^if) having one double and one simple zeros 

Gonsider the graph Figure 5.(n). Here /i is a simple zero and /2 is a double zero. The real 
solution occurs between fi and /2. At /i we have P^ifi) = > 0, so graph of the 

function / is concave up at .^i. At double zero /2, / —t /2 as ±cxd. Hence we have a 

solitary wave solution with amplitude fi — f 2 <0. The other case gives unbounded solution. 


To sum up, we can give the following proposition for £ = 3 case. 


Proposition 2.1 Equation (1.9) may admit solitary wave solutions when the polynomial 
function P^^f) admits (i) one double and three simple zeros (ii) two double and one simple 
zeros (Hi) one triple and two simple zeros (iv) one quadruple and one simple zeros (v) one 
double and one simple zeros. 


3 Exact Solutions by Using the Chebyshev’s Theorem 


The Ghebyshev’s Theorem is given as following |10j . 












Theorem 3.1 Let a, b, c, a, [3 be given real numbers and aj3 Q. The antiderivative 

I = [ x°‘{a + (dx^Ydx (3.1) 


is expressible by means of the elementary functions only in the three cases: 

(1) + c e z, (2) e z, (3) c e z. 


(3.2) 


The term x°'{a + is called a differential binomial. Note that the differential binomial 

may be expressed in terms of the incomplete beta fnnction and the hypergeometric function. 
Let us dehne u = /3x^/a. Then we have 

1 2+i+co-^T5 /1 + fl \ 

I = -a » » Byi^—,c-iy 

1 g+l I „ g+l l+g ^ f O ~\~ 3 ^ iTuT^ 


-a 


i+a^/a + l ^ i + a + o \ 
[3 0 u b - - - ;uj, 


_ (3.3) 

1 + a Vo 0 / 

where By is the incomplete beta function and F{t, k] rp, u) is the hypergeometric function. 


Our aim is to transform the system (1.1) to {y') = Pi+ 2 {.y) by taking f{Cj = ^ + ay[Yf) 


We can apply Chebyshev’s Theorem to this equation if we assume that P^^ 2 {y) reduces 
to the form ^^+2(2/) = Ay~‘^°‘{a + (3yY~‘^'^, where —2a — 2c + b = i + 2. For i = 3, let 


u{x,t) = /(O = 7 + (y.y{l3f) in (1.9), then the equation becomes 

(y'f = P-,(y) = aiy^ + + ayp + ayy + Oo, 


(3.4) 


where 


a\ — 

02 = 

03 = 

04 = 

«5 = 

Og = 


a^ 

1 / 5«^7 , 0-2 V 

^(_ + 3ac). 

+ 6 dc^ — 2adi + 12007 ^, 

1 


-p ^ - Icdi + Ad 2 + 18 c 7 ^ - 6^17 + 67^ + 4 c^ + 180^7^, 


1 


(^Sds + 12c7^ H—^ + 18 c^ 7 ^ — Scdi'y + 80^7 + 80(27 ~ 60(17^ j, 

5 

_2^2' (^^ 31 ' + + Y + 80(4 - Acdi ' f '^ - 2c(i7^ + 4 c^ 7 ^ + 4 c( 27 ^) • 


a /32 

1 


To apply the Chebyshev’s Theorem 3.1 we assume that P^iy) reduces to the following form, 

Ps{y) = Cy-^^ {a + fiyT^Y (3.5) 


where a, b, c are the constants in Theorem 3.1, b — 2a — 2c = 5, and C is a constant 


Here we present the cases mentioned in the Proposition 2.1 Other cases are given in 
Appendix A. 
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1) Let P^{y) = y{a + This form corresponds to the case of one simple or one simple 

and two double zeros. We have y' = ■P.y^P{a + (3y‘^) so 


y 


+ (3y‘^) ^dy = ±.^ + A. 


(3.6) 


Here a = —1/2, 6 = 2, and c = —1. Hence 

(1) ^^ + c=-3/4^Z, (2) ^^ = 1/4^Z, (3) c=-leZ. 


For (3), from (3.6), by letting u = (dy'^/a we obtain 

jti) =±^ + A 


(3.7) 


Here choose a = —1, /d = 1, the integration constant A = 2, plus sign in (3.6), we have 

arctanh(y^) + arctan(y^) = —+ 2), (3.8) 

and the graph of this solution for ^ < —2 is given in Figure 6 



Figure 6: Graph of a solution of {y'Y = 2/(—1 + y'^Y 


This is a kink-type solution. 


2) Let P^y) = 2/(*^ + fdyY- This form corresponds to the case of one simple and one 
quadruple zeros. We have y' = ■iiy^P{a -1- jdyY so 


y 


^/^(a -h (5y) ‘^dy = -|- H. 


(3.9) 


Here a = —1/2, 6=1, and c = —2. Hence 

(1) -1-c =-3/2 ^ Z, (2) ^^ = 1/2^Z, (3) c=-2gZ. 


For (3), from (3.9), by letting u = (dy'^/a we obtain 

■1 




Vv 


a{jdy -f a) 


arctan ( 

+-= + a (3.10) 


a 
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Here choose a = 1, /? = 1, the integration constant A = 2, plus sign in (3.9) we have 


l + y 


+ arctan(y^) = ^ + 2, 


(3.11) 


and the graph of this solution for —2 < < vr/2 — 2 is given in Figure 7 



Figure 7: Graph of a solution of (y')^ = y(l + y)'^ 


3) Let Psiy) = y^(a + /3y^). This form corresponds to the case of one simple and one 
double zeros. We have y' = ±?/(a + l3y'^YP so 

j y~^{a + /3yY~^^^dy = + A. (3.12) 

Here a = —1, 6 = 3, and c = —1/2. Hence 

(1) ^^ + c=-l/2^Z, (2) (3) c=-l/2^Z. 


For (2), from (3.12), we obtain 


a/ ,o/3 ..\ \\V3 


y = (g(‘““*'^ (^^“('^±0) - 1 


(3.13) 


Here choose a = 4, /3 = 4, the integration constant H = 2, and plus sign in (3.12) we 
obtain 

l/(0 = -(sech(3e + 6))2/3 (3.14) 

and the graph of this solution is given in Figure 8 



Figure 8: Graph of a solution of {y'Y = ?/^(4 + 4?/^) 
This is clearly a solitary wave solution. 
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4) Let P^iy) = y‘^{a + Py)^. This form corresponds to the case of one double and one 
triple zeros. We have y' = ±y{a + (3y)^P so 

j y~\a + l3y)-^/^dy = + A. (3.15) 

Here a = —1, 5=1, and c = —3/2. Hence 

(1) + c = -3/2 ^ Z, (2) = 0 e Z, (3) c = -3/2 ^ Z. 


For (2), from (3.15), we obtain 

2 


ayj (dy + a 


2arctanh { ASjP 


Here choose a = 1, /? = 1, H = 2, and plus sign in (3.15) we have 

. — 2arctanh(A/?/ + 1) = ^ + 2, 

V2/ + 1 

and the graph of this solution is given in Figure 9 



(3.16) 


(3.17) 


Figure 9: Graph of a solution of {y'Y = y‘^{l + yY 


5) Let P^iy) = y'^{a + (3y). This form corresponds to the case of one simple and one 
quadruple zeros. We have y' = Py'^{a + /3yY^'^ so 


-2 


y 


(a + M + A. 


(3.18) 


Here a = —2, 5=1, and c = —1/2. Hence 


, , n -)- 1 , , ^ 

(1) ——+ c=-3/2^Z, 


For (2), from (3.18), by letting u 


-a ‘f( - 1, ^;0;ii) = 


(2) ^ = -l6Z, (3) c = 

Yy^!a we obtain 

/^arctanh(^^) _ 
ay 


-1/2 ^ Z. 


±e + A (3.19) 
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Figure 10: Graph of a solution of {y'y = 2/^(1 + y) 


Take a = 1, /3 = 1, A = 2, and plus sign in (3.18) we get 


arctanh('^j/ + 1) — = .^ + 2, 


y 

and the graph of this solution for ^ > —2 is given in Figure 10 


(3.20) 


4 An Alternative Method to Solve (/')^ = Pn{f)^ When 
n > 5: Factorization of Polynomials 

When we have (/')^ = Pn{f) for the reduced equation it becomes quite difficult to solve 

such equations for n > 5. For this purpose, we shall introduce a new method which is based 

on the factorization of the polynomial Pn{f) = P£+2{f), n > 5. Let the polynomial P£+2{f ) 

£+2 


have i + 2 real roots i.e. P£+ 2 {f) = B n(/ — fi), B is a. constant. Define a new function p{^) 

i=l 

so that / = /(p(0)- Hence we have {ff = . By taking 


Ni 


'0 - 


2=1 

N2 




2=1 


where iVi + -^2 = ^ + 2,i? = up, we get a system of ordinary differential equations. Solving 
this system gives the solution of the degenerate coupled f'-KdV equation. 

For illustration, we start with a differential equation where we know the solution. Con¬ 
sider 

{y'y = a(y - yi){y - y2){y - y?i){y - 2 / 4 ). ( 4 . 1 ) 

Let y = y{p{0) so 

^ a{y - yi){y - y2){y - y?i){y - yi)- 


(dy^. 

£dp\ 

\dp} ' 
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Take 


J) = Av - V\)(v - Vi) 
= v(y - vz){y - Vi), 


(4.2) 

(4.3) 


where k/x = a. We start with solving the equation (4.2). We have 

dy 


Viv - yi)iy - y^) 


dy 


b^ -1 


= ±y/Kdp, 


where a = {yi + y2)/2 and b = {yi—y2)/2. Let {y — a)/b = cosh.^. This gives us dz = ±y/kdp. 
After taking the integral we obtain 

arccosh^^^-^—j + Ci, Ci constant, 

y = h{p) = a + b cosh.{^/Kp + C 2 ), C 2 constant. (4.4) 


so that 


Now we insert this solution into Eq.(4.3). Let cosh.{^/Kp + C 2 ) = v. Hence we get 

= + C 3 , 


dv 


i/nby/v"^ — — A)^ — B'^ 


(4.5) 


where A = ( 2 / 3 + 1 / 4 -|/i-|/ 2 )/( 2 /i- 2 / 2 ), S = ( 1 / 3 -|/ 4 )/(|/i- 1 / 2 ), and C 3 is a constant. Solving 
(4.5) and using (4.4[) gives 


^ //i(//3 - y2) + y2{yi - i/3)sn^((i/2)^/t(i/3 - //2)(//4 - yl){^JJ^i + C'a), k) 

^ (?/3 - 2 / 2 ) + {yi - |/3)sn2((l/2)^K(j/3 - j/ 2)(?/4 “ ?/i)(V+^ + C's), k) 

which was obtained in [8]. Here sn is the Jacobi elliptic function and k is the elliptic modulus 
satisfying P = {f^ - - / 3)/(/2 - /3)(/i - A). 

Now we apply our method to the case when the polynomial Pn{f) is of degree n = 5. We 


have several possible cases, but here we shall give the case when the polynomial (1.9) has 
five real roots. Other cases are presented in Appendix B. 


Case 1. If (1.9) has hve real roots we can write it in the form 

iff = a{f - /i)(/ - fi){f - h){f - h){f - h), 


(4.7) 


where /i, / 2 , fs, / 4 , and /s are the zeros of the polynomial function P 5 {f). Now define a new 


function p{^) so that / = f{p{^)). Hence (4.7) becomes 




(4.8) 
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i. Take 


= -(f - fiKf - h){f - h){f - fi) 

tf = -“(/-A)- 


(4.9) 

(4.10) 


In [H] we have found the solutions of Eq.(4.9). One of the solutions is 

. , /.(/a - /i) + fzih - /4)sn"((l/2) V(/i - h){ f 2 - h)p, k) 

(A - /i) + (/i - /4)sn^((l/2)v^(/i - h)Ui - fi)p, *;) 

where k is the elliptic modulus satisfying = (/a — / 2)(/4 — /i )/(/3 — /i )(/4 — / 2 )- We use 


(4.11) 


this solution in the equation (4.10) and get 


dp 


A + Bsn‘^{up, k) 


>0 Vf5 - h{p) Jo \l C + Dsn^{up,k) 


dp = + Cl, Cl constant, (4.12) 


where 


A = fs-fi,B = fi-U,C= (/s - mfs - fi), D = {h- h){fi - U), 
and w = (l/2)y^(/i - - Ia)- For particular values; /^ = /a = -1, /g = 3, /4 = 0, 


/s = 2, a = 1, Cl = 0, and choosing plus sign in (4.12) we get the graph of the solution f{^) 
given in Figure 11 



Figure 11: Graph of the solution (4.11) with (4.12) for particular parameters 


This is a solitary wave solution. 
ii.Take 

ijp) ^ - fi)if - h)if - fo) 

(^) ^ h(/-/4)(/-/5), 


(4.13) 

(4.14) 
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where = a. Eq.(4.13) has a solution 


/ = Kp) = (/2 - /i)sn^((l/2)\//3 - + C'l), k) + /i, (4.15) 

where k satisfies = (/i — / 2 )/(/i — /s) and C'l is a constant. We use this solution in 


Eq.(4.14) and get 

/■^ dp 


dp 


= ±A//i^ + C2, 
(4.16) 


Jo a/( h(p) - fi){h{p) - /s) Jo Asn.^{up + 5, /c) + Bsn^{u:p + 5, /c) + C 
where C 2 is a constant and 

-4 = (/2 - /l)", B = (/2 - /l)(2/l - /4 - /s). C = (/i - /4)(/, - /s), = (1/2 )Jk(/3-/i). 

and (5 = (l/2)V/3 - /iCi. For particular values; fi = 1, /2 = 2, /s = 3, /4 = -1, /s = -2, 
K = /i = 1, Cl = 0, and choosing plus sign in (4.16) we get the graph of the solution f{^) 
given in Figure 12 



Figure 12: Graph of the solution (4.15) with (4.16) for particular parameters 


This solution is periodic. 


iii. Take 


|)' = '^(/-/l)(/-/2) 

= M/- AX/- AX/-/s). 


where np = a. Consider first the equation (4.17). 


(4.17) 

(4.18) 


We have 


df 


V(f - SOU - A) 


= -^y/kdp. 


(4.19) 
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This equality can be reduced to 


df 


= ±y/Kdp, 


(4.20) 


b^l(^] -1 


where a = (/i + f 2 )/‘^, b = {fi — /2)/2. By making the change of variables (/ — a)/b = u, we 
get 


du 


J \/u^ — 1 

where Ci is a constant. Thus we obtain 


— i-\/Kp + Cl, 


f - ^ ^ 

arccoshu = arccosh( — ) = ±y/ Kp + Ci, 


(4.21) 


(4.22) 


which yields 

f = h{p) = a + bcosh{^/Kp +Cl), a = (/i +/2)/2, b={fi-f2)/2. (4.23) 


Now we use this result in (4.18|), 
rp 


dp 


'o \/{.h{p) - h){h{p) - f4){h{p) - h) 

dp 


A cosh^(^/Kp + Cl) + B cosh^(y^p + Ci) + C cosh(A/K/) + Ci) + D 
= + ^* 2 , ( 4 . 24 ) 

where C 2 is a constant and 

A = b^, B = &^(3a-/3-/4-/5), C = b{{a-f 4 ) + {a-f^) + {a-fi){a-f^)}, 


and D = {a- / 3 )(a - / 4 )(a - /s). 

For particular values; /i = 5, /2 = 1, fz = fi = 


-3, /s — —4, K — 4, /i — 1, Cl — C 2 — 0, 


and choosing plus sign in (4.24), we get the graph of the solution f{^) given in Figure 13 



Figure 13: Graph of the solution (4.23) for particular parameters 
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iv. Take 


1)^ - 

'dp^ 2 


(^) ^ Kf - f2){f - f3){f - f4){f - f5), 


(4.25) 

(4.26) 


where up = a. Consider the Eq.(4.25). We have 

df 




= ±v^dp. 


(4.27) 


Integrating both sides gives 

/ = h{p) = (^ ± + Cl j + fi, Cl constant. 


(4.28) 


We nse this resnlt in the eqnation (4.26) 

r dp 


>0 V W) - h){h{p) - h){h{p) - h)ih{p) - /s) 

dp 


^(f ± ^Cip + 2l)(f p2 ± ^Cip + 5)(f p2 ± ^Cip + C)(fp2 ± ^Cip + D) 

= ±a/7^^ + 6*2, (4.29) 

where C 2 is a constant and 

A = Cl + fi- /2, B = Cl + h- /s, C = + /i - /4, D = C' + /i - /s. 

For particnlar valnes, /i = 1, /2 = /s = 2, /4 = 3, /s = 4, k = 4, /i = 1, and Ci = C 2 = 0, 


and choosing pins sign in (4.29) we get the graph of the solution f[^) given in Figure 14 


This is a solitary wave solution. 



Figure 14: Graph of the solution (4.28) for particular parameters 
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5 Conclusion 


We study the degenerate ^-coupled KdV equations. We reduce these equations into an ODE 
of the form (/')^ = P£+ 2 (/), where Pf +2 is a polynomial function of / of degree i + 2, i > 3. 
We give a general approach to solve the degenerate Pcoupled equations by introducing 
two new methods that one of them uses Chebyshev’s Theorem and the other one is an 
alternative method, based on factorization of P£+ 2 (/), i >3. Particularly, for the degenerate 
three-coupled KdV equations we obtain solitary-wave, kink-type, periodic, or unbounded 
solutions by using these methods. 
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APPENDIX A. Other Cases Using Chebyshev’s Theorem for Degenerate Three- 
Coupled KdV Equation 

Here we present other cases that we use Cheyshev’s theorem to solve degenerate three- 
coupled KdV equation. 

1) Let P^iv) = y{a -I- This form corresponds to the case of one simple zero or three 
simple zeros. We have y' = ■P.y^P{a -|- (3y^YP so 

j y~^'‘^{ci + (3y^)~^/‘^dy = + A. 

Here a = —1/2, 6 = 4, and c = —1/2. So we have 

(1) + c =-3/8 ^ Z, (2)^^ = 1/8^Z, (3)c=-l/2^Z. 

b b 

Hence we cannot obtain a solution through the Chebyshev’s theorem. 

2) Let Pdiy) = + Py"^)- This form corresponds to the case of one triple or two simple 

and one triple zeros. We have y' = ±y^P{a -|- (dy'^yP so 

j + yy^-^/'^dy = + A. 

Here a = —3/2, 6 = 2, and c = —1/2. So we have 

(1) ^^ + c=-3/4^Z, (2) ^^ = -l/4^Z, (3) c=-l/2^Z. 

6 6 

Hence we cannot obtain a solution through the Chebyshev’s theorem. 
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3) Let P^iy) = y^{a + Py)"^- This form corresponds to the case of one double and one 
triple zeros. We have y' = ±y^P{a + (3y) so 


y + M 'dp = + A. 


(6.1) 


Here a = —3/2, 6=1, and c = —1. Hence 

(1) ^^ + c=-3/2^Z, (2) ^^ = -l/2^Z, (3) c=-leZ. 


For (3), from (6.1), by letting u = fiy^/a we obtain 

1 „ 1 


/ I I V 9 2^/5 arctan f 

- + ^' (8-2) 


^y/y 


4) Let P^iy) = (a + I3y)^. This form corresponds to the case of a zero with multiplicity 
hve. We have y' = ±{a + Py^P so 


(a + f3y) ^Pdy = + A. 


(6.3) 


Here a = 0, 6 = 1, and c = —5/2. Hence 

(1) + c = -3/2 ^ Z, (2) = 1 e Z, (3) c = -5/2 ^ Z. 


For (2), from (6.3), by letting u = /3y^/a we obtain 

9 _ \ 2 


a ^uf(i, 2;2;m) = 


Here we get 


y = 


3/3(/]y + a)T2 


— a 


= ±.^ + A. 


(6.4) 


( 6 . 5 ) 


/?n9/?2(±J + 5l)2, 

5) Let P^iy) = (a + /3y^). This form corresponds to the case of one simple zero. We have 
y' = ±{a + fiy^yP so 

Ji^ + I3y^)~^^^dy = ±^ + A. 

Here a = 0, 6 = 5, and c = —1/2. We have 

(1) + c =-3/10 ^ Z, (2) ^^^ = 1/5^Z, (3) c=-l/2^Z. 

0 b 

Hence we cannot obtain a solution through the Chebyshev’s theorem. 


APPENDIX B. Other Cases Using Alternative Method for Degenerate Three- 
Coupled KdV Equation 
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Here we present the other cases for the method of factorization of the polynomial P 5 {f)- 


Case 2. If (1.9) has three real roots we can write it in the form 

(/')" = a(f - /i)(/ - Mf - A)(of + bf + c), 


( 6 . 6 ) 


where /i, / 2 , and /a are the zeros of the polynomial function P^if) of degree hve and 


6^ — 4ac < 0. Now define a new function p(^) so that / = /(p(0)- Hence (6.6) becomes 

Wr- = (f )'(f)' = “(/ - /i)(/ - f)(/ - A)(“/' + V + c), (6.7) 


i. Take 


(?) ^ - fi)if - f2){f - h) 


•.dp 

(?) ^ /^(«/^ + ^/+ c), 

where np = a. Consider the hrst equation above. We have 

, = Py/ndp. 

V(f - /i)(/ - k)(f - A) 

We know from Case l.ii that Eq.( |6.8 ) has a solution 

/ = ^(p) = (/2 “ /i)sn^((l/2)\//a — fiiy/np + Ci)) + /i, 


( 6 . 8 ) 

(6.9) 


( 6 . 10 ) 


( 6 . 11 ) 


where Ci is a constant. We use this solution in (6.9) and we get 
rp 


dp 


dp 


— Py/p^pC2^ ( 6 . 12 ) 


Jo yjah^^p) + hh{p) + c Jo ■\/o4sn"^(a;p + 5) + Hsn2(a;p + 5) + C 
where is a constant and 

^ = o (/2 — /i)^, -S = (2a/i + &)(/2 — /i), C = {af"^ + bfi + c), u: = (1/2)^(/s ~ /i), 

and 5 = (l/2)/^^C'i. 


ii. Take 


Pf = f/-/i)(/-A) 

= p(/ ~/3)(«/^ + ^/+ c), 


(6.13) 

(6.14) 


where np = a. Consider first the equation (6.13). It has a solution 


f = h{p) = ai + bicosh{y/Kp +Cl), m = (/i +/2)/2, bi = {fi - f2)/2, 


(6.15) 
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as it is found in Case l.iii. We use this solution in (6.14) and get 


dp 


'0 - /3)(ah2(p) + bh{p) + c) 


dp 


A cosh^(^/Kp + Cl) + B cosh^(^/Kp + Ci) + C cosh(^/Kp + Ci) + D 

= + C 2 , (6.16) 

where C2 is a constant and 

A = abl, B = bl{ 3 aai + b — /s), C = bi{ 3 aal — 2001/3 + 2 &Oi — f^b + c), 
and D = {ai — /3)(o^o + bai + c). 


iii. Take 


dfY (f 

'dp'^ 2 


(^) ^ - f2){f - f3){af+ bf+ c), 


where np = a. Consider the equation ( |6.17 ). It has a solution 

/ ^/i(p) ^ ( ± ^p + Cl) +/i, Cl constant 


(6.17) 

(6.18) 

(6.19) 


as it is obtained in Case l.iv. We use this result in Eq.(6.18) 

r _ dp _^ 

Jo V W) - h)W) - h){ah^{p) + bh{p) + c) 

r dp 


sj (f± ^Cip + p4)(f/)2 ± y^c'ip + ± <^c^p3 + cp^ + Dp + E) 

= + ^2, (6.20) 

where C 2 is a constant and 

A = Cl + h- / 2 , B = Cl + h- /3, C = K{b + QaCl + 2o/i)/4, 

and 

D = ±^Ci{2af, + 2aCl + b),E = {h + Cl){ah + aCl + b) + c. 


iv. Take 

(^) ^ + ^/ + c) (6.21) 

(^) ^ Kf - h){f - f2){f - fs), (6.22) 
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where Kfi = a. Consider the hrst equation above. We have 


df 


/±s 

w 


b \ 2 


— + Cl, 


+ 1 


where W = /‘la. Let jW = tan 6 '. Hence the above equality becomes 


sec 6d6 = y/a{^y/kp + Ci), 


which gives 


In I sec^ + tan6^1 = ^/a{l!l^/Kp + Ci). 

Hence after some simplihcations and taking the integration constant zero we get 

/ = h{p) = ±A/4ac — smh{^/aKp)/la. 


We use this result in (6.22) 
rp 


dp 


Jo ^/W) - fi){h{p) - f2)(h{p) - fs) 

r _ dp _ 

J^ \JA sin\^{^/akp) + B sinh^(A/aKp) + C sinh(A/aKp) + D 
= + C2, 

where C 2 is a constant and 

A = ±(4ac- b^)~^^^/8a^,B = (b^ - 4ac)(fi 4-/2 + / 3 )/ 4 a^ 

and 

C = ±V4ac-b‘^{fif3 + 72/3 + /i/2)/2a, D = -/i/2/3- 


(6.23) 


(6.24) 


Case 3. If (1.9) has just one real root we can write it in the form 

(/O^ = ~ /i)(/^ + + <22/^ + aif + oo), 


(6.25) 


where fi is the zero of the polynomial function P^if) of degree five and the constants a*, 
i = 0 , 1 , 2 ,3 are so that + 03 /^ + 02 /^ + aif + oq 7 ^ 0 for real /. 


Now define a new function p(^) so that / = /(p(0)- Hence (6.25) becomes 

(/O^ = (df) ^ 4“ “ 3 /^ + 02 /^ + cq/+ flo)- (6.26) 
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i. Take 


dfY ff 

Tp) = 


+ <^ 2 /^ + ^if + < 20)5 


where k^ = a. Consider the hrst equation above. It has a solution 

a/k n 2 


(6.27) 

(6.28) 


/ = h{p) = ( ± + <^ 1 ) + /i, Cl constant 


as it is obtained in Case l.iv. We use this result in (6.28) and get 
rp 


dp 


'0 ^/h^{p) + a^h^ip) + a2h‘^{p) + cliHp) + ao 


= + C 2 , 


(6.29) 


(6.30) 


where C 2 is a constant. Here we do not get a simpler expression than the original equation. 
Hence it is meaningless to use the method of factorization of the polynomial for this case. 


ii. Take 


dp) 

d^Y 

d^) 


— ^(/^ + ^ 3 /^ + ^ 2 /^ + fli/ + flo) 
= p(/-/i), 


(6.31) 

(6.32) 


where up = a. Here for simplicity, we will use the form /'^ + 6/^ + c instead of the polynomial 
function of degree four above. Since f^ + bf‘^ + c = (/^ + 6/2)^ —6^/4 + c, we will take c > 6^/4 
to get an irreducible polynomial. If c = &^/4, then take 6 > 0. From ) =«:(/4 + 6/2 + c) 


df 


= ±a/kp + Ci, Cl constant. 


0 .//4 + 5J2 + C 


We have 


f{p) = 




\/—h + a/— 4c + 


sn(( a/2/2) (Cl ± V^p) \J-h + a/-4c P), 


where k = ^ ^ — • Then we use this function / in (|6.32|) and get 


\fAdp 


"'0 y/A/2csn( (Ha/2/2) (Cl ± i/^p), k) - Afi 
where A = \/—b + a/— 4c + P. 


= ±^/p^ + C 2 , C 2 constant, (6.33) 
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Now take c = 6^/4 with 6 > 0. In this case we have 


df 


/o P + b/2 


= ±y/Kp + Ci, Cl constant, 


which gives the solution 


/ = Pp) = i^tan (^^(V^p + Ci)]. 


(6.34) 


Now use the equation (6.32). We have 
rp 


dp 


° ■/±^tan ( ^(a/kP +Cl)) -/i 


where C 2 is a constant. Solving the above equation gives 


— + ^* 2 , 


±4 fv^.2, \{Ai- A^f + A\ , 

- . 9 , . { ——AAn — -——--j -Oarctan 

pJib{Al + AI)AP 4 3 l{Ai + A2y + All 


2^1^3 


Al-Af + Al 


I — =tA/h'C + C' 2 , 


where 


Ai = ]J ± 2 y^tan {^^{^/Kp + Cl)) - 4/1, 2I2 = \j 


+ 26 - 2 / 1 , 


and A^ — ij 4/^ + 26 + 2^. 
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